DSF 21/2002 
astro-ph/0301438 



Neutrino energy loss rate in a stellar plasma 



(N 



S. Esposito, G. Mangano, G. Miele, I. Picardi, and O. Pisanti 

O 

Dipartimento di Scienze Fisiche, Universitd di Napoli "Federico II" 

and 

Istituto Nazionale di Fisica Nucleare, Sezione di Napoli 
Complesso Universitario di Monte S. Angelo, Via Cinthia, 1-80126 Napoli, Italy 
E-mail: sesposito@na.infn.it, mangano@na.infn.it, miele@na.infn.it, picardi@na.infn.it, 

pisanti@na. infn.it 

> ' 
OO 

m 

O ! Abstract 

m 

We review the purely leptonic neutrino emission processes, contributing to the 

■ energy loss rate of the stellar plasma. We perform a complete analysis up to the 
first order in the electromagnetic coupling constant. In particular the radiative elec- 
tromagnetic corrections, at order a, to the process e + e~ — > vv at finite density and 

> ■ temperature have been computed. This process gives one of the main contributions 

to the cooling of stellar interior in the late stages of star evolution. As a result of 
the analysis we find that the corrections affect the energy loss rate, computed at 

■ tree level, by a factor (—4 4- 1)% in the temperature and density region where the 
pair annihilation is the most efficient cooling mechanism. 



PACS number (s): 13.40. Ks, 95.30. Cq, ll.10.Wx 



1 Introduction 



One of the crucial parameters which strongly affect the stellar evolution is the cooling 
rate. Stars during their life can emit energy in the form of electromagnetic or gravitational 
waves, and/or as a flux of neutrinos. However, in late stages a star mainly looses energy 
through neutrinos, and this is pretty independent of the mass of the star. In fact, white 
dwarfs and Supernovae, which are the end points for stars with very different masses, have 
both cooling rates largely dominated by neutrino production. An accurate determination 
of neutrino emission rates is therefore mandatory in order to perform a careful study of 
the final branches of star evolutionary tracks. In particular, a change in the cooling rates 
at the very last stages of massive star evolution could sensibly affect the evolutionary 
time scale and the iron core configuration at the onset of the Supernova explosion, whose 
triggering mechanism is still lacking a full theoretical understanding pQ . 

The energy loss rate due to neutrino emission (hereafter denoted by Q) receives con- 
tribution from both weak nuclear reactions and purely leptonic processes. However for 
the rather large values of density and temperature which characterize the final stages of 
stellar evolution, the latter are largely dominant. The leading leptonic processes are the 
following: 

i) pair annihilation e + + — > v + V 

ii) z/-photoproduction 7 + e ± — > e ± + v + V 

iii) plasmon decay 7* — > v + V 

iv) bremsstrahlung on nuclei e ± + Z — > e ± + Z + u + 17 

Each process above results to be the dominant contribution to Q in different regions of 
the core density-temperature plane. For very large core temperatures, T > 10 9 °K, and 
relatively low density, p < 10 5 g cm -3 , the pair annihilation is the most efficient cooling 
process. For the same values of densities but lower temperatures, 10 8 °K < T < 10 9 °K, 
the z/-photoproduction gives the leading contribution. These density-temperature ranges 
are the typical ones for very massive stars in their late evolution. Finally, plasmon decay 
and bremsstrahlung on nuclei are mostly important for large (p > 10 6 g cm -3 ) and 



extremely large (p > 10 9 g cm -3 ) core densities, respectively, and temperatures of the 
order of 10 8 °K <T < 10 10 °K. Such conditions are typically realized in white dwarfs. 

Starting from the first calculations of Ref.s j21 El > a systematic study of the energy loss 
rates for processes i)— iv) has been performed in a long series of papers [2]- [22]- in all these 
analyses the pair production rate i) has been evaluated at order G%, i.e. at the zero-order 
in the electromagnetic coupling constant a expansion, whereas the remaining processes 
ii)-iv) are at least of order aG\. Thus to correctly compare the energy loss rates for all 
processes i)— iv) it is worth computing QED radiative corrections to pair annihilation rate 
i), which may lead to a sensible change in the cooling rate Q. This is the aim of our 
analysis, whose results have been briefly reported in Ref. [22] • in this paper we give all 
the details of the calculations. 

The paper is organized as follows. Section 2 is devoted to a brief summary of Born 
amplitude calculation for pair process, while in Section 3 we report the details of order a 
QED corrections. Neutrino photoproduction and plasmon decay are discussed in Section 
4 and 5, respectively. Our results are summarized in Section 6. 



2 Born amplitude for pair annihilation process 



a = e, p, t, and the 4-momenta are defined aspi^ = (£a,2, Pi, 2) and = (^1,2, ^1,2)- The 
energy loss rate induced by this process is obtained by integrating the squared modulus 
of the invariant amplitude M e + e -^ UaVa over the phase-space of the involved particles, and 
summing over the flavour of final neutrinos, 



for e with temperature T and degeneracy parameter £ e , and in J2 sp in.a a sum over an 
particle polarizations and final flavours is performed. Notice that, as long as neutrino 
mean free path is large enough that they can leave the star without any further inter- 



Let us consider the annihilation process e (pi) + e + (p2) 



v a (qi) + M<?2), where 




.(2.1) 



The quantities F±(E) = exp j|J±£ e j + 1 



are the Fermi-Dirac distribution functions 



2 




action, there is no relevant neutrino component in the stellar plasma. Thus no neutrino 
distribution function is present in the expression for the energy loss rate. 

In this Section we evaluate Q e + e - in the Born approximation (hereafter denoted with 
Qf +e -), i.e. in the limit of a four-fermion electroweak interaction and no electromagnetic 
radiative correction (see Figure^). Let us consider first the pair annihilation in electron 
neutrinos. At first order in perturbation theory, the diagram in Figure contributes to 
the invariant amplitude M e + e -^ VeVe with two terms coming from W ± and Z° boson ex- 
change, respectively. By using the low energy expression for the vector boson propagators, 
one has 1 

MY +e -^ aVe = -_^- 5 [«( ft ) 7M (l - l5 )u( Pl )] [v(p 2 mi - 7i)t;(fe)] , (2.2) 

and 

M? +e -^ UaVa = - f - 2 [v(P2h,(C v - C Al5 )u( Pl )} [rZ(gi)7 /i (l-75Mg 2 )] , (2.3) 

O COS £/\y ''vz 

1 We use natural units, Ti = c = k = 1. 
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where C v = 2 sin 2 6 W - 1/2, C A = -1/2 and a = e, //, r. By using a Fierz transformation 
on Eq. ()2.2|) and summing it to ()2.3|) one gets the total amplitude in Born approximation, 

MX-^v, ~ TsMfc)] HP2mC v ' - C A ' l5 )u( Pl )} , (2.4) 

with 

C v ' = 1 + C V = \ + 2sm 2 6 w , (2.5) 

CV = 1 + C A = ~ . (2.6) 



The squared modulus of the amplitude (|2.4j) . summed on the polarizations of the incoming 
and outgoing particles, can be expressed as the product of the two tensors, Tf£) and Tf$, 



EI^WJ 2 = f ^r^r . (2 . 7) 



M", _ _ p = — — i ^ i ^" 

e+e — >v e v e I Q 

spin 

The tensors T^) and Tj$ in Ea.(l2~Tl) can be both decomposed in a symmetric (S) and 
antisymmetric (A) part in the indices fiv, namely 

rp(e,u) _|_ rp(e,u) 

S £u) s g + ^ ( (2 _ g) 



T 

In their product only the SS and AA combinations survive, but the latter disappears after 
performing the integration over neutrino phase spaces. This integration can be performed 
by using the Lenard formula, namely 

IP^f ^ {4) (P ~qx- V = ^(2pV + rt 2 ) e(p°) &(p 2 ) , (2.10) 
J 2 ui\ J 2 uj 2 24 

where denotes a generic 4- momentum. By means of Eq. (j2.1U|) . the quantity in curly 
brackets of Eq. ()2.1j) . but only for electron neutrinos, takes the form 



( B + - I 



1 fd^rd\± p2 _ gi _ 

2tt) 2 J 2wi J 2u 2 ^ n 
^ K + Pl ■ p 2 ) [(C v 2 + C2) Pl -p 2 + (2 C{ 2 - C' 2 ) m 2 ] . (2.11) 



For Up or v T production, the W exchange term is instead absent and thus only the neutral 
current contributes. In this case, by using the previous arguments the same expression 



is obtained but with the substitution C V ',C A ' — > C V ,C A . By virtue of the above 
results, the total Qf +e -, obtained by summing on neutrino flavour, reads 

lo7T° JO E\ JO -C/2 

c , 2 UE\E\ | 9E 1 E 2 E\ + El | \ | c;2 /4gg iff + g 

V \ TXl\ W? e m e / \ m e m e 



(2.12) 



where we have denoted with Cy A = (l + CV ij4 ) 2 + 2C 2 A and we have performed the angular 
integrations. Note that Qf+ e - depends on the temperature T and the electron degeneracy 
parameter £ e only. 

It is customary to recast the dependence of the energy loss rate on £ e (or the electron 
chemical potential) in terms of the matter density, p, the temperature, T, and the electron 
molecular weight, p e , 

^ = £*if' (2-13) 

Pe i s±i 

where in the above expression the sum is performed over all nuclides, Zi and Ai stand 
for the atomic and the weight number of the 2-nuclide, respectively, and X.% is its mass 
fraction. To this aim, by requiring the electrical neutrality of the plasma, we have 

n e - -n e+ = N A — , (2.14) 

Pe 

where n e ± are the number densities and Na is the Avogadro number. The degeneracy 
parameter £ e can be then obtained by inverting Eq. (j2.14|) . namely 

vr^kll. Bf^™ U { - _ f. } + x - „p { . 1 1. } + J • (215) 

Once performed numerically the integration over electron/positron momenta in Eq. (|2.12j) . 
the energy loss rate due to pair annihilation in the Born approximation can be expressed 
as a function of T and pj p e only. The results are shown in FigureEl where Qf+ e - is plotted 
as a function of pj p e , for the following values of temperature T = 10 8 , 10 8 ' 5 , 10 9 , 10 10 °K. 

3 Radiative corrections to pair annihilation 

A consistent computation of the aG 2 F corrections to the energy loss rate induced by pair 
annihilation is obtained by considering in addition to the tree level graph of Figure EJi, 
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Figure 2: The energy loss rate due to pair annihilation process in the Born approximation, 
Qf+ e -, is here plotted for T = 10 8 , 10 8 - 5 , 10 9 , 10 10 °K (from bottom to top). 




the radiative diagrams of Figures [Tj3-[TJi. The corrections of order aG 2 F are then obtained 
via their interference. 

It is worth while observing that since the typical energy carried by the outgoing neu- 
trino pair is at most of the order of 1 MeV, one can safely neglect the electroweak radia- 
tive corrections to the four-fermion effective interaction (involving additional weak boson 
propagators), and only consider the gauge-invariant set of purely QED contributions. 

When the processes we consider take place in stellar interiors, we must take into 
account the influence of the electromagnetic plasma (namely, e ± ,7) when computing 
radiative corrections. To this end we have performed calculations by using the Real Time 
Formalism for finite temperature quantum field theory. In this framework the thermal 
propagators for electrons and photons read as follow 



iS(p) 
iD aP (k) 



{j> + m) 



p 2 — ml + ie 



- r F (p) 



~9af3 



k 2 + ie 



(3.16) 
(3.17) 



with 



T F (p) = 2tt 5 (p 2 - ml) [F_(bo|)e(p ) + F + (\p \)Q(-p )] 
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(3.18) 



T B (k) = 2nS(k 2 )B(\k \) 



(3.19) 



where 0(x) is the step function and B(x) is the Bose-Einstein distribution function. The 
first terms in Eqs. (j3.16|) and (|3.17|) are the usual T = Feynman propagators, while those 
depending on the temperature (and density) through the distribution functions describe 
the interactions with real particles of the thermal bath. 

According to the different diagrams contributing to radiative corrections, it is custom- 
ary to classify these corrections as follows: 

- electron mass and wavefunction renormalization (Figures H^-^); 

- electromagnetic vertex correction (Figure [Hi); 

- 7 emission/absorption (Figures HMD 1 )- 

3.1 Zero— temperature radiative correction 

By using in the evaluation of radiative corrections of Figures ^d-^I the first term only 
in the propagators of Eqs . (|3 . 1 . (|3. 17jl . one obtains the off-shell contribution to the ra- 
diative amplitude, and the corresponding result is known as zero-temperature or vacuum 
correction. 

This issue has been addressed in Ref. [21] for the process e~ + v — > e~ + u, but we 
can easily obtain the desired corrections for pair annihilation by using crossing symmetry. 
The vacuum correction to be inserted into Eq. (j2.1j) is then 

spin, a 

where 



7T 



aG F m e u 2 5(E 2 ,lu 2 ) 



(3.20) 



5(E 2 ,oo 2 ) = g 2 L iv 1 (E 2 )+V 2 (E 2 



z-l + 

+ fH \ \\(E 2 )(\ :T f V 2 (E 2 ) 

'm P z 



m e z 



2u 2 



2-1 + 



m e z 



+ 9 LR UVi(E 2 )-V 2 (E 2 ) 



u 2 



2u 2 
-2V 2 (E 2 



z-\-z A 



, (3.21) 



V 1 (E 2 )=Re<;(2 1ogY 



En 



2|p 2 



log 



E 2 + 


|F>2 1 


E 2 — 


|F>2 1 



E^ 

|F»2 I 



Li 2 



|P2 




E 2 


- m e 


2 


P2 
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£ 2 H 


- |F»2 1 


E 2 - 


- |F»2 1 


sin 2 


9w + 



W^H-^logli^T) , (3-23) 

with 

g\ = 3sm A 6 w -sm 2 9 w + ^ , (3.24) 

g 2 R = 3sin 4 # w , (3.25) 

g LR = sin 2 9 W (3 sin 2 9 W - ^ , (3.26) 

z = (E2 + m e )/uj2 and A is a small photon mass introduced to regularize the infrared 
divergences. Remind that the dilogarithm function Li 2 (x) is defined by 

Li 2 (x)^-/W g(1 - t) . (3.27) 
Jo t 

The expression ()3.2U|) depends on the infrared regulator A and diverges for A — > 0. This 
is quite obvious since as it is well known, when computing QED radiative corrections, 
in addition to the pair annihilation process one must consider also the bremsstrahlung 
radiation accompanying the process (Figures ^-0i). In fact, only the combination of 
virtual and real photon corrections is free from infrared divergencies. 

The bremsstrahlung process will be considered in detail in Section 13.2.31 but here we 
anticipate the strategy for a careful and reliable numerical computation of the integrals 
in Eq. (j3.33J) . In fact, numerical cancellation of divergencies is very hard to handle and 
we prefer to divide the bremsstrahlung contribution into a "soft" part (according to the 
photon energy being lower than some fixed threshold e) which will cancel the infrared 
divergencies of the pair process, and a "hard" one which, in our case, will depend on the 
thermal photon distribution function. The last one will be considered in Section 13.2.31 
while we now report the Soft-Bremsstrahlung (S.B.) squared modulus. Again we use the 
result of Ref. [21] with a suitable crossing transformation and thus we get 

E \AM^ VaV f = ^aGWe^I^e) \gl + g\ (1 - zf + g LR (™) 



spm,a 



(3.28) 
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where 



- L 





= Re j 


E 2 


- IP2I 


E 2 


+ |F»2 1 



E% . / E 2 + Ip2 
log 



2 p 



log 



and 



E 2 + 


P2I 


E 2 — 


P2I 


L(x) 





1 - 2 log 



E 2 - IP2L 
^ P2I 



+ 



Ei 



2|Pa| 
+ 1 - 2 log 2 



£2 + 


|Pa| 


£ 2 - 


IP2I 



^log]l^ 

t 



(3.29) 



(3.30) 



Remarkably, by summing the expressions (|3.2(J|) and (|3.28|) . since the term [Vi(E 2 ) + 
I-y(E 2 , e)] does not depend on the infrared regulator A, the total squared amplitude is now 
infrared divergence free. Then, the term in curly brackets in Eq. (J2.1j) . once integrating 
over all variables but one by using the 5-function, becomes: 
1 r d 3 qx f d 3 q 2 



$ = 



(2tt) 



2lu, 



2uo 



■(5 w (pi + 



P2-qi- Q2) 



E (|AMJ + = £ °__ q77 J 2 + |AM ( 



e+e _ — >!/ a ^ a I 



/2 



16 aGp-mf r^M uj 9 
dUJ2 \^\ 

n -IWrn IP2I 



x \5(E' 2 ,u' 2 ) + I,(E' 2: e) 



2,2/1 A2 , / 

V W 2 



(3.31) 



where we denote with a prime the quantities in the electron rest frame, and the integration 
limits are 

E' 2 + m e E' 2 + m e 



uj m = m e 



E 2 + |p' 2 | + m e ' ~ " E' 2 - \p' 2 \ + m e 

Using this result, the correction to the energy loss rate at zero-temperature is 



(3.32) 



T=0 

e+e — 



32 vr 4 Jo 



d\pi\ d\p 2 



d(cos9 12 ) | Pl | 2 |p 2 | 2 E l + ^ 2 ^F-(Ei) F + {E 2 ) 



Ei E 2 



(3.33) 



where 9i 2 is the angle between pi and P2, $ stands for $ boosted to the comoving frame 
and the integrals must be numerically evaluated. 

3.2 Thermal radiative corrections 

The true thermal radiative corrections to the Born estimate of the neutrino pair pro- 
duction process come from considering in the evaluation of radiative diagrams of Figures 
^i-EJl the terms involving at least one thermal part of the propagators (J3.16|) and (J3.17)) . 



These contributions do not involve ultraviolet divergencies due to the presence of the Fermi 
and Bose function in the thermal propagators. Nevertheless, the results corresponding 
to each of these corrections are not free from infrared divergencies and, in general, are 
gauge-dependent too. The inclusion in our calculations of the photon emission/absorption 
diagrams in Figures ^-Qi is then required to overcome these difficulties, as we have di- 
rectly and carefully checked. In the intermediate steps, however, we have to regularize 
the actual divergences and, to this end, we have explicitly subtracted all divergent terms 
expanding the squared amplitudes in a Laurent series around the pole singularities. This 
method follows quite closely what has been already used in Refs. [2E]- 

3.2.1 Mass and wavefunction renormalization at finite temperature 

As already stated, at order a, the thermal mass shift and thermal wavefunction renor- 
malization corrections come from the interference of the diagrams in Figures with 
the tree-level one of Figure [IJi, having subtracted the zero-temperature contribution. 
The thermal mass correction to neutrino energy loss rate may be obtained by replacing 
mass m e with the renormalized value m^± = m e + 5m± in the Born expression for the 
rate and subtracting the zero-temperature limit 1)2.12)1 . The self-energy for an electron of 
energy E and momentum p has been calculated in Ref. [26 a and we refer the reader to 
this paper for details. At order a the thermal mass shift for e ± is given by 

an m0 a r°° „, , Ik 



air a /"°° k 

6m± = T 2 + / rf|k| -J- (F±(E k ) + F T (E k )) 

Sm e nm e Jo E k 

am e 

2tt|p 



am f°° Ik I 

+ d|k|U (F ± (E k )\ogC„ + F T (E k )logC + ) , (3.34) 

In p Jo E k 



m e 2 + 


Pi 


|k| ± EE k 


m e 2 — 


IpI 


|k|± EE k 



C± = % " ' , — - , (3.35) 



with 



where E k = J |k| 2 + m 2 and E = y |p| 2 + m 2 . By using the thermal renormalized masses 
in the expressions for energies appearing in the Born rate in Eq. (j2.12j) . after some algebra 
we obtain the following thermal mass correction to the neutrino energy loss rate 

AQX- = AQ££ + AQ^ e F + + AQ^I , (3.36) 
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with 



AQf + f_ 



aG 2 F 
9tt 6 



rf| Pl |rf|p 2 |rf|k|| Pl | 2 |p 2 | 2 |k| (^r + ^r 



x F_(E 1 )F f (^)B(|k|)^A + A) 



(3.37) 



X 



18^6 7 



d|pi|d|p 2 |d|k| 

,2 / 



|Pi| 2 |P2| 2 |k| 



Ek 



l- + l.)F_(t\)E + [E 2 )E ± (E, ) 

j±i\ Hi 2 , 



|k| \E X + E 2 ) + 2 \^i|pi 



J> logCii + ^^logC^ 



-E2IP2I 



(3.38) 



m P 



{C' A 2 + Of) 3 + 6 



,E\E 2 



+ 



3(C. 



: 2 - 2c; 2 ) 



(cj? + c; 2 ) 

E\E 2 \ 



3(Cj 



' £i#2 / 
3^1 + 



E\E 2 



Eo 



m e 



|Pl| 2 |P2| 2 

mi 



1 



(3.39) 



/a 



{c? + c; 2 ) 3 + 6 



+ 



£1 

l-F+(^ 2 ) 



,E\E 2 



nit 



- 3(C; 2 - 2C, 



— m P 



Si 



E 2 {Ei + £"2) 



3(C 



: 2 - 2c; 2 ) 



E\E 2 \ 



E\ En ( E-\ En 

3^^ [1 + -^ 



+ 



mi 



nit 



|Pi| 2 |P2| 2 

mi 



1 



mi 



(3.40) 



where C±i j2 are obtained from C± in Eq. (j3.35|) with the substitution (E, p) — > (Ei <2 , p± <2 ). 
Note that in Eq. (J3.36|) we have subtracted all divergent terms as described above. 

Diagrams in Figures ^>-^: are also responsible for the thermal wavefunction renor- 
malization correction which, in the calculation for the neutrino energy loss rate, can be 
obtained by using a thermal renormalized projector on positive/negative energy states 
as described in Refs. (213 EEj- The projector on positive energy states for an electron of 
4-momentum p = (E, p) can be cast in the following form 

p 1 + mf + 5 + 



E 



2E 
(jj) + m e )A 



iPl 



m P 



B 



(3.41) 
(3.42) 



with p = (0, p) and where the momentum- dependent functions A and B are given by 



A 



a 



2nE 



d|k||k| 



2S(|k|) 



log 



E + 


IpI 


E- 


IpI 
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+ E_(E k ) 



2|k| 



XE-E k 



+ 



logC_ 



- F + (E k ) 



2|k| 



IpI 



B 



a 



2tt|p| 3 £ 



\(E + E k y 

f°° Ikl / 

/ d|k|y £j4|p|£-2m 2 log 
Jo E k I \ 



E + 


IpI 


E — 


IpI 



(3.43) 
S(|k|) + (2|k||p|E 



- m 2 e (E - E k ) logC_) F_(E k ) + (2|k||p|£ - m 2 e (E + E k ) \ogC + ) F+(E k )] , (3.44) 



The projector on negative energy states reads instead 

j> — mf + <5_ 



A 



R 



2E 

(j)-m e )A + \ f + 



iPl 



m e 



B 



(3.45) 
(3.46) 



and A and B can be obtained from A and £> by replacing £ e — > — £ e . The contribution 
due to e thermal wavefunction renormalization is then obtained by using the thermal 
renormalized projectors A^ in the evaluation of the Born rate. This procedure, having 
subtracted the contribution due to the mass renormalization, gives at order a 



AQ 



G 2 F 

127T 5 



jTdlPildlPal |pi| 2 | P2 | 2 (J- + Jr) F_{E X )F + {E 2 



(C 2 + CI 2 ) (E 1 E 2 + mfj {2E,E 2 + mfj + - | Pl | 2 |p 2 | 



(Ai + A 2 ) 



(3.47) 



+ 2 - (C 2 + C^ 2 ) | Pl | 2 |p 2 | 2 (£i + B 2 ) - 2 (2 Ci 2 - C't) (E X E 2 + ml) (Ipil 2 ^ + |p 2 | 2 £ 2 )) 

The functions Ai t2 , B± j2 , A± j2 , B\ )2 are obtained from A, B, A, B with the substitution 
(E,p) —* (Ex j2 ,pi j2 ). Again, in the numerical computation of the integrals appearing in 
the expression for AQ^+ e - we have subtracted all divergent terms. 

3.2.2 Vertex renormalization at finite temperature 

The lowest order electromagnetic vertex correction to the neutrino energy loss is pro- 
vided by the interference term between the diagram in Figure and the tree amplitude 
of Figure [T^i. In the amplitude of the vertex diagram, three particle propagators ap- 
pear, each of them consisting of two terms according to Eqs. (|3.16|) and (J3.17|) . However, 
the computation of the thermal vertex renormalization correction to the rate is simpli- 
fied by 4-momentum conservation arguments. In fact, the mentioned three propagators 
would produce eight terms in the amplitude, one of which is a pure vacuum term already 
considered in Sect. 13.11 Three of the remaining terms, proportional to a Fermi-Dirac 
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distribution function times a Bose-Einstein distribution function, give no contribution to 
the rate since they correspond to an electromagnetic vertex with a real photon and two 
real massive particles, which is not allowed by 4-momentum conservation. Moreover, the 
term proportional to two Fermi functions does not contribute to the interference term 
with the tree level amplitude since it is purely imaginary. Therefore, we are left with only 
three terms, one proportional to a Bose-Einstein function and the others proportional to 
a Fermi-Dirac distribution 

A *- = -Sgj* I ****** + ±) F^)FAE 2 ) 

r B (k) r F {k + Pl ) r F {k- P2 ) 



X 



[(k + pi) 2 — m 2 ] {(k — p 2 ) 2 — m 2 ] [(k — P2) 2 — ml]k 2 [(k + pi) 2 — rn^k 2 



(3.48) 



where 



$ y = 512 Cy 



m 2 



+ Pi - Pi) (-2 Pi ■ P2 (2m 2 e +p l -p 2 ) - k-p 2 (3m 2 e + 2p 1 ■ p 2 
+ k ■ pi (3m 2 e + 2 k ■ p 2 + 2 p l ■ p 2 



- 512 m\{k-p x ) 2 + m 2 e (k-p 2 ) 2 + 2k-p 2 pi-p 2 (m 2 e + p ± ■ p 2 ^ 

- 2k-p 1 p l -p 2 {m 2 e + k-p 2 +p 1 -p 2 ^+2p 1 -p 2 [~m A e + (p l -p 2 ) 2 ^j . (3.49) 

The infrared-safe part of these terms (in the limit A; — > 0) is obtained by expanding the 
matrix element function <£>y in powers of k. 

For the Bose-Einstein term (the first one in square brackets in Eq. (j3.48|0 . due to the 
presence of the function B(\ko\), one has to subtract from $y the terms in powers of k of 
order and 1. By choosing the frame where the momentum k lies along the z-axis and 
Pi is in the x-z plane, denoting with 8\ y2 the angle between k and pi 2 and with <fi the 
azimuthal angle of the vector p 2 , after some algebra one gets 

x F_(E 1 )F + ( J E 2 ) J B(|k|)/ B , (3.50) 

where 0:1,2 = cos 6*1,2 and Ib is reported in Appendix 1X1 

The fermionic part (the second and third terms in square brackets in Eq. (j3.48j) ) con- 
tains, instead, the expression 

6(k 2 + 2k. Pl ) 

dk (k 2 -2k- P2 )k 2 ® v ' (3 - 51) 
13 



which, assuming that the integration in dk Q does not introduce divergent terms, is pro- 
portional to 

d 3 k- 



v 



k ■ (pi + p%) k ■ pi 
By using the properties of the 5-function, for k — > we can write 



(3.52) 



-E x + J E% + |k| 2 + 2|k|| Pl |a;i 



iPiFl 



E\ - |pi| 2 xi 



E x 1 1 2Ef 1 1 

and substituting in the expression in (|3.52|) we find that 



O Ik 



d A k 



5{k 2 + 2k- Pl ] 



d\k\ 



(3.53) 



(3.54) 



(k 2 -2k-p 2 )k 2 ~" |k| 
Thus, the infrared-safe part of the fermionic term is obtained by subtracting from $y 
the terms in powers of k of order only. With the same notations as above, after some 
algebra we get 

*V,F 



6144^ Jo d \*wP2m £ dx xdX2 |Pl1 ^ |k| + ^ 



x F^{E l )F + {E 2 ) [I Fl F^(E k ) + J F2 F+(£ fc )] 



(3.55) 



where E^ - 

/ ^ 

10 



'Fx 



k| 2 + m 2 . The functions 2 are defined as follows 

+ 



\{k — pi — p 2 ) 2 - m 2 ] (k - pt)' 



k =E, 



[(k -px- p 2 ) 2 - m 2 ](£; - p x ) 



ip 2 



2- 



[{k + px + p 2 ) 2 - m 2 ] (k + p 2 ) : 



+ T 

k =E k \- 



+ Pi + P2) 2 - m 2 e ](k + p 2 y 



ko=~E k 



k =-E k 



(3.56) 



(3.57) 



where &y ) is obtained from $y subtracting the terms of order in k and making 
the substitution k — > k — px (k — > k + p 2 ). The integrations in Eqs. (|3.56|) . (|3.57|) may be 
performed analytically, since they involve integrals of the form 



2tt 



nx + n 2 cos (ft + n 3 cos 2 (ft + n 4 cos 3 ■ 



(3.58) 



10 dx — d 2 cos (f> 

where n^, dj do not depend on 0. However, since the denominator above may vanish 
for given values of the momenta involved, the result of the integrations depends on the 
kinematical region one considers. For each function I Fl and Ip 2 we distinguish three 
different cases: 
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the integrand functions are just polynomials in cos(0) and no pole is involved 

7r(2n!+n 3 ) 



di 



(3.59) 



2. d 2 ^ 0, 



> 1: 

inside the integration region in no pole may occur and the result is 

72 



-2vr 



U2 

d 2 2d 2 



n 4 d\ n 3 d{ ri 4 



d\ 



d\ 



l(h + d 2 1 
d\ — d 2 d\ + d 2 



ni H ; 1 s — h 



(3.60) 



3. d 2 ^ 0, 



< 1: 



in this case the integrand function may develop a pole and the integral I has to be 
evaluated in the principal value sense: 

-7r [2 d\ d 2 n 3 + 2 d\ n 4 + d\ (2 n 2 + n 4 )] 



d3 



(3.6i; 



The final expressions of the functions Ip x 2 are reported in Appendix [X] Thus, the 
thermal vertex renormalization correction to neutrino energy loss rate due to pair anni- 
hilations reads 

AQ e V = AQ?£_ + AQX- . (3.62) 



3.2.3 Bremsstrahlung: 7 emission/absorption 

As stated above, in order to eliminate the infrared divergencies present in the radiative 
diagrams of Figures ^i-EJl it is necessary to include the rates of processes where a photon 
is either absorbed or emitted (see Figures [T^-Eft). The energy loss rates due to 7 emission 
(absorption) are given by the sum of squared amplitudes of the processes of 

Figures HMD 1 ; namely 

= p/^/^/S (Bl+£2 - |k|)F - (£l)F+(B2)[1 + B(|k| 

■5 A (p 1 +p 2 -qi 



x 



I /"Ah. f d 3 q 2 xi; 



(2 



2wi 



2w 2 



-q 2 -k) Yl \M, 

spin,c 



e+e-—>i/ a i/ a j\ 



(3.63) 
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AQ 



A 

e+e" 



X 



^ V y ± ^ spin, a 



1 efe _ 7-n'a!/ Q 



(3.64) 

Note that in the computation of the amplitudes for both processes, only the vacuum part 
of the electron or positron propagator should enter since, otherwise, we would have an 
electromagnetic vertex with a real photon and two real electrons or positrons which is not 
allowed by 4-momentum conservation. 

The term proportional to unity in the statistical factor (l+5(|k|)) entering in Eq. (j3.63J) 
is responsible for the spontaneous emission in vacuum. For small |k| this contribution has 
been already discussed in Section \'d.l\ fSoft-Bremsstrahlung). hence we have to consider 
only the hard-photon emission part corresponding to photon energies above the introduced 
cutoff e. 

For the 7 emission, the term in curly brackets in Eq. ()3.63|) is given by 



aG 2 F 
24 



Q(E 1 + E 2 - |k|) Q(m 2 e + Pl ■ p 2 - k ■ (p! + pa)) 



(pi • k) 2 

while for 7 absorption we have 



+ 



v 12 



(p 2 -k) 2 (pi-k)(p 2 -k) 



(3.65) 



24 



+ 



+ 



$f 2 



(3.66) 



_{pi-k) 2 (p 2 -k) 2 {pi - k)(p 2 - k)_ 

The functions $f , $f , are reported in Appendix El while for the photon 

absorption are obtained from the corresponding <3>f by simply replacing the 4-momentum 
k with —k. The B-functions appearing in Eq. (j3.65|) . coming from the Lenard formula 
used in performing neutrino momentum integration, give the kinematical conditions for 
the process to occur. Obviously, for the 7 absorption we have no kinematical constraints 
and we have omitted the superfluous B-functions 2 . 

In the limit |k| — > the integrand function in Eq. (j3.63|) (where the expression in 

Eq. fl3.65j) has been substituted) shows some divergences which cancel the corresponding 

2 In this case the arguments of the 0-functions are always positive for all values of the particle 4- 
momcnta, as it can be easily checked, for example, in the electron-positron center of mass frame. 
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infrared singularities in the radiative corrections considered previously. To pick up the 
divergent terms, we note that the relevant expression appearing in AQf +e - is the following: 

H (E 1 + E 2 -\k\) <1> E [1 + B(\k\)] . (3.67) 

Thus, by expanding <& E in powers of |k|, the infrared-safe part is obtained by subtracting 
the terms up to order in |k| to the expression proportional to the Bose function. As 
far as the factor which does not contain -B(|k|) is concerned, we observe that no change 
has to be performed because we must integrate this term for |k| ranging from the cutoff 
e to the infinity, since the soft part of the (vacuum) bremsstrahlung has been already 
considered in Section I3~T1 

Analogous divergent terms for |k| — * appear in the integrand function for the energy 
loss induced by photon absorption. In this case the relevant expression to be considered 
is 



and, as in the previous case, the infrared-safe part is obtained by subtracting terms up to 
order in |k|. 

Finally, the integrations in the electron, positron and photon 3- momenta in Eq. (j3.63J) 
must be performed by taking into account the 0-functions in Eq. (j3.65j) . We choose the 
frame where the photon momentum k lies along the z-axis and pi is in the x-z plane. We 
also denote with 9i j2 the angle between k and pi j2 and with the azimuthal angle of the 
vector p 2 . 

The kinematical constraints can be implemented as follows. The integration range for 
the variables |pi|, | p» 2 1 is not limited (they run over the entire [0, oof interval), as well 
as those for 9\, 9 2 (ranging from to 7r). The integration in the modulus of the photon 
momentum |k| is bound to the region |k| < E\ + E 2 , thus implementing the condition 
coming from the first B-function in Eq. 1)3.65)) . The remaining kinematical constraint can 
be written as: 



k| 2 



(Ei + ^ + lkl) $ A B(\k\) , 



(3.68) 



Oi COS0 < CL 2 , 



(3.69) 



with (s i)2 = sin 9 lj2 , x 1>2 = cos#i i2 ) 



a i = IP1IIP2IS1S2 



(3.70) 
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a) b) 

Figure 3: Feynman diagrams for the neutrino photoproduction processes. 



a 2 = E X E 2 - |pi||p 2 ki^2 - |k|(^i + E 2 - |pi|xi - |p 2 |x 2 ) + m 2 e . (3.71) 

This gives no effective condition when a± = 0, a 2 > or a\ ^ 0, a 2 /ai > 1, thus leaving 
the integration over <fi ranging from to 2tt, while Eq. fl3.69|) cannot be fulfilled when 
a± = 0, a 2 < or ai ^ 0, a 2 /ai < —1, thus resulting in AQ^+ e - = for these cases. 
In the remaining cases, that is a\ ^ 0, — 1 < a 2 /a\ < 1, the integration interval in is 
[a, 27r — a] with a = arccos(a 2 /ai). 

For the photon absorption case, since there are no kinematical constraints, the inte- 
gration on |k| as well as on the angle is not restricted to a given region but ranges over 
all possible values. 

The final expression for the sum of the neutrino energy loss rate induced by photon 
emission and absorption is the following: 

AQf +e _ + AQf +e - = 35^7 / d\p 1 \d\P2\d\Mdx 1 dx 2 l ^^ F^(E 1 )F + (E 2 ) 

x[I e0 Q(E 1 + E 2 -\k\) + (I e e(E 1 + E 2 -\k\) + I a )B(\k\)} , (3.72) 

where the function J e o, I e , I a are reported in Appendix iBl 

4 Neutrino photoproduction 

Neutrino photoproduction processes j(k) + e ± (pi) — > e ± {p 2 ) + v a (qi) + v a (q 2 ) contribute 
to the neutrino energy loss rate, at the same perturbative order as the electromagnetic 
corrections to the pair annihilation process. Thus, they must be included in a comprehen- 
sive study of the star cooling rate. The relevant diagrams for photoproduction on e~ or 
e + , shown in Figure El are simply obtained by crossing the corresponding ones for photon 
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absorption in Figures^, [TJi, and the neutrino energy loss rates take the following forms 



l 7e — — >e~ v a v a \ 



(4.73) 

/ / 4? ^ + ^ - Pi - 9l - 92) X! \ M ie+^e+u a V a 

1 2 spin, a 



(2vr) 2 



(4.74) 



Differently from what happens for the corrections to the pair process, analyzed in the 
previous Section, in the calculations of photoproduction rates we must take into account 
the photon effective mass m 7 depending on the density and temperature of the plasma. 
This effect could be neglected in all the e + e~ — > vv radiative corrections because it 
yielded an irrelevant change in the neutrino energy loss rate, but it is important for the 
photoproduction processes which exhibit a strong dependence on the background density 
in the high density region. 

In the following we consider a massive photon with the following renormalized mass 

E3 

m\ = — / d|k| JJ- (F + (E k ) + F_(E k )) , (4.75) 

' 7T JO 

where we have adopted the same notations of Eq. (j3.34J) and the completeness relation 

E« A) = -V + ^. (4-76) 

Note that a finite photon mass eliminates infrared divergencies due to the Bose distribution 
function in the rates. For the sake of brevity, in the following we describe the calculations 
for photoproduction on electron only. The corresponding ones for photoproduction on 
positron can be obtained by replacing £ e — > — £ e . 

The term in curly brackets in Eq. (j4.73|) is given by 
8 a G 2 

%h = -f ®{E l -E 2 + E 1 )Q(2m 2 e + m 2 y -2p l -p 2 + 2k-{p l -p 2 )) 

]>'' _|_ P" _|_ Ph rjrj^ 



(m 2 + 2pi • A;) 2 (m 2 — 2p 2 -k) 2 (m 2 + 2pi • /c)(m 2 — 2p 2 ■ k) 
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The squared amplitude function $p fe , §^ h , $ p |, are reported in Appendix O The 9- 
functions appearing in Eq. ([4.77)1 . coming from the Lenard formula used in performing 
neutrino momentum integration, give the kinematical conditions for the process to occur. 
The remaining integrations on the incoming electron and photon and outgoing electron 
3- momenta in Eq. ()4.73j) are performed, in close analogy to what done in Section I3.2.3| 
by choosing the frame where the photon momentum k lies along the z-axis and pi is in 
the x-z plane (denoting with # 1)2 the angle between k and pi 2 and with the azimuthal 
angle of the vector p 2 ). In particular the integration in |pi|, |k|, 6\, # 2 ranges over all 
possible values for these variables, while the integration field for | p»2 1 is limited by the first 
6-function in Ea. (l4~7Tl) : 

< |p 2 | < y /\p 1 \ 2 + E* + 2E^E 1 =p^ ax ■ (4.78) 

The remaining kinematical constraint, which has to be used to perform the integration 
on the angle 0, can be written as: 

61 cos0 > b 2 , (4.79) 

with (s 1>2 = sin^ lj2 , x 12 = cos 6^2) 

61 = 2|pi| |p 2 | sis 2 , (4.80) 
b 2 = 2E l E 2 -2| Pl | \p 2 \x 1 x 2 -2E 7 (E 1 - E 2 ) + 

+2 |k| (| Pl | x x - |p 2 | x 2 ) -Iml-rn^ . (4.81) 

Note that the condition K7fy is not fulfilled when bi = 0, b 2 > or bi ^ 0, 62/61 > 1, 
thus resulting in Q ie - = 0, while it is trivially realized for 61 = 0, 6 2 < 0, where case 
the integration field for is [0, 27r]. The constraint (|4.79|) can be non trivially satisfied 
only for 61 7^ 0, — 1 < 6 2 /6i < 1, when the integration region for is [0, 0\ U [2n — (3,2tt] 
with j3 = arccos(6 2 /6i). The final expression for the energy loss rate induced by neutrino 
photoproduction on electron is the following 

aG\ too rrf- * |pi| 2 | P2 | 2 

Qie- = Qn70 7 / o pi ok / a p2 / — 

J072 7T' JO JO J-l ki\ki2 

x F^E^l-F-iE^BiEjI^ , (4.82) 
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V 

Figure 4: Feynman diagram for the plasmon decay into a neutrino- antineutrino pair. 

where the function I p h is reported in Appendix O 

The corresponding expression for the rate induced by neutrino photoproduction on 
positron is easily obtained from Eq. (j4.82|) with simple substitutions, 

aG 2 F f°° [P™ ax f 1 I F»i 1 2 1 F>2 1 2 

Qie+ = Qn70 7 / d|p2|d|k| / d|pi| / dx x dx 2 — 
dU7/7r' jo Jo J-i kt\ht2 

x F + (E 2 ) [1 - F+iEk)] B(Ej) 4 , (4.83) 

and I' h is obtained from I ph with pi <->• p 2 3 . Thus by using the expressions ()4.82|) and 
fl4.83jl . one gets the total energy loss rate due to neutrino photoproduction as 

Qje Q~/e~ ~t~ Q~fe+ • (4.84) 

5 Plasmon decay 

The plasma process 7*(/c) — > v a {qi) + V a {q 2 ) is of higher order in a with respect to the pair 
annihilation process (it is of order otyfaG 2 F ) but, nevertheless, its contribution to neutrino 
energy loss rate is dominant for a wide region of stellar temperatures and densities [21. 
Then, in this Section, we discuss this process and evaluate the corresponding energy loss 
rate. 

The decay 7* — > v a + V a takes place only if the photon momentum is a time-like 4- 
vector such that > |k|. To see when this constraint applies, we have first to calculate 
the dispersion relation for a photon in a thermal plasma and this is recalled in Appendix 
ID1 For transverse (T) propagation modes there is no effective constraint, while for lon- 
gitudinal (L) ones the decay occurs only for |k| < |k| max , where |k| max is the maximum 

3 Note that the kinematics of the process is the same as for the photoproduction on electron by replacing 
the 4- vector p\ with p2- 
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plasmon momentum in the bath whose expression is reported in Appendix \D\ 

At lowest order the process 7* — > v a + V a is described by the Feynman diagram in 
Fig. E]and the loss rate is given by 

Q r = Q L r + Q T r , (5.85) 

where the partial rates 

vLT 1 r d 3 k 



>< (5 ' 86) 
correspond to longitudinal and transverse photon mode propagation, respectively. The 
computation of the invariant amplitude for the considered decay can proceed through an 
effective photon-neutrino interaction 

M^ uV = u( qi )-f(l - 75)^2) T a ,(k) e a (k) , (5.87) 

where e"(k) is the polarization 4- vector for longitudinal (A = L) or transverse (A = T) 
thermal photons and, by considering only the thermal on-shell propagation of e ± in the 
loop in Figure |H Thus one gets 

eGV f d 4 p 



T F {p-k) | r F (p) 



(5. 



San = (P~ k )aPv + ip~ k )^Pa + P • kg afl , (5.89) 

A afl = e^Wtf . (5.90) 

The quantity in curly brackets in Eq. (j5.86j ) is thus given by [19J: 

^ T = - ^I'^ ^P^ , (5-91) 

PT = -Z L Q a0 V a "(v^Y , (5.92) 

Pf = -2Z T R aP Y a ^ (V pv X , (5.93) 



R aP = g aP - - Q al3 , (5.95) 
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Figure 5: The total radiative corrections normalized to the Born approximation result for 
the pair annihilation process for T = 10 8 , 10 8,5 , 10 9 , 10 10 °K (from top to bottom). 



where u a = (1,0) is the medium 4- velocity and we have used the orthogonality condition 
K a T aiM = 0, with K = (Ej, k) and E 1 = \J\k\ 2 + rri*. Note that for longitudinal and 
transverse photon modes in a plasma we have [2E1 Ell 

E e A e? = - Q aP Z L , (5.96) 

A=0 

E e " el* = -2K#Z T , (5.97) 

A=±l 

respectively, where the k— dependent functions Z^t can be found in Appendix |Dj The 
energy loss rates are, then [T5J |2T] 

Or = ^^/ lk|maXrf l k ll k l 2fi (^L)^n| , (5.98) 

Q T r = d\k\\k\ 2 B(E jT )Z T W T , (5.99) 

with n L T reported in Appendix iDl 

6 Results and conclusions 

In this paper we have presented an exhaustive computation of the energy loss rates in 
neutrinos of the stellar interior. A consistent analysis of all leptonic processes up to order 
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a in the electromagnetic fine structure constant has been performed: neutrino energy 
loss rate due to pair annihilations, evaluated in Born approximation, Qf+ e - ()2.12j) . has 
been corrected including both vacuum and thermal radiative corrections, the latter being 
computed in the real time formalism. By using Eqs.flZZE), (jHHJ), flUSHl), flUHI, flH3E|), 
(|H.fi2|) . and (|H.72jl . the total radiative correction AQ e + e - results 

AQ e + e - = Qe+e- ~ Qe+e~ 

= aqt + =°- + aqX- + Ag^ e - + aq £ V + Agf +e . + aqX- ■ (6.100) 

In Figure we plot the ratio AQ e + e -/Q^ +e ^ as functions of the plasma density for 
some values of the temperature, namely T = 10 8 , 10 8 ' 5 , 10 9 , 10 10 °K. The corrections are 
found to be of the order of few percent and negative for high temperatures, implying 
that for these temperatures the energy loss is sensibly decreased. For fixed temperature, 
AQ e + e - /Qf+ e - goes to a constant value for low density. This can be easily understood, 
since in this limit the plasma is weakly degenerate, and therefore the energy loss rate 
depends on temperature only. At large densities the ratio AQ e + e -/Qf +e _ decreases and 
reaches larger negative values. However, for such high densities the pair annihilation rates 
are exceedingly small and thus this process gives only a marginal contribution to the star 
cooling. 

In order to perform a comprehensive study, we have also recalculated the contribution 
from i/-photoproduction and plasmon decay processes. The corresponding numerical re- 
sults are presented in Figure where we show the several contributions to total neutrino 
energy loss rate. When possible, we also show for comparison the results of Ref. |2"2"] . 
In particular we find a fair agreement with the results reported in literature. For com- 
pleteness we also show the energy loss rate due to bremsstrahlung on nuclei, denoted with 
Q e z, which have been produced using the analytic fitting formula of Ref. |22j . 

In summary, Figure El shows that, as well known, for large temperatures and not too 
high densities, pair annihilation dominates over the other three processes, while for low 
densities z/-photoproduction dominates over plasmon decay and bremsstrahlung on nuclei. 
On the other hand, for large densities the most relevant process is plasmon decay, whose 
rate however, along with those of all other processes, rapidly falls down for extremely 
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Figure 6: The energy loss rate versus p/fx e due to pair annihilation including radiative 
corrections Q e + e - , photoproduction Q ie and plasmon decay Q 7 * (solid lines) for several 
temperatures (see for definitions Eqs. (|2.1|) . (|4.84|) . (|5.85j) ). The dotted lines refer to the 
analogous results of Ref. [22], where the rate for bremsstrahlung on nuclei is also com- 
puted. The effect of AQ e + e - to pair annihilation cannot be appreciated in this logarithmic 
scale but for the largest temperature 10 10 °K. 
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Figure 7: The regions in the T — p/fi e plane where each of the processes i)-iv) contribute 
for more than 90% to the total energy loss rate. We also show the contours for the relative 
correction AQ/Q% (see text) for the values 1%, 0%, -1%, -2%, -3%, -4%. 

high densities. This is a genuine plasma effect as noted in [23] • Consider, for example, 
the behaviour of z/-photoproduction energy loss. As already noted in [Hj, the decrease for 
very large densities is achieved only if one consistently takes into account the increasingly 
large photon thermal mass. In fact with a massless photon the z/-photoproduction curves 
in Figure |H1 would rather reach a constant value. The main effect of is a lowering of 
the values of the Bose distribution function for photons, i.e. a smaller number of thermal 
photons. This reduces the energy loss rate induced by z/-photoproduction. 

In Figure[3we show the regions in the temperature-density plane where a given process 
contributes to the total energy loss rate Qxot = Qe+e- + Q^e + Q 7 * + Qez (including 
radiative corrections to pair annihilation) for more than 90%. We also summarize there 
our results on the radiative corrections to pair annihilation processes, by plotting the 
contours corresponding to AQ e + e - IQ\ ot = 1%, 0%, —1%, —2%, —3%, —4%, where Q% t ls 
the total emission rate with pair annihilation calculated in Born approximation. These 
contours lie almost entirely in the region where the pair annihilation process gives the 
main contribution to the total energy loss rate. This result may affect the late stages of 
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evolution of very massive stars by changing their configuration at the onset of Supernova 
explosion. 
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A Vertex correction functions 



With the same notations of Section 13.2.21 the function I B appearing in Eq. (|3.5U|) is given 
by: 

1 r 2n 



'b 



A Jo 



&v 8 



pi ■ kp 2 ■ k 



+ 



fc =|k| 



Pi ■ kp 2 ■ k 



feo= 



256C; 2 vr [2E 1 E 2 -2\p 1 \\p 2 \x 1 x 2 -m: 



E x - 


|Pi 


Xi 


E 2 — 


|P2 


x 2 



E 2 - 


|P2 


x 2 




|Pl 


Xi 



+512 CI 2 n (m^ + S lJ B 2 -|p 1 ||p 2 | 



X\ x 2 



(A.l) 



where $y 5 is $y once subtracted the terms of order and 1 in k. 



On the other side, after performing the integration involved in Eq. (J3.58|) . one obtains 
in the three cases of Section 13.2.21 



1. d 2 = 0; 



If 



x 
x 



-128 7T 



(m? e — k- pi) [k- (pi +p 2 ) - EiE 2 -m? e + |pi| |p 2 | x 1 x 2 ] 



{c; 2 [| P i| 2 |p 2 | 2 (-1+a;! 2 ) (-l+x\ 
(Q Ex E 2 - 2 k ■ p 2 + 5 m 2 - 6 |pi| |p 2 | 



x x x 2 



+ 2 [Ex E 2 + m 2 e - \px\ |p 2 | xx x 2 ) (2 E{ E 2 -5k-p 2 m 
+ k ■ pi (2 k ■ p 2 + 3 m 2 ) - 3 - 3 m 2 |pi| |p 2 | x x x 2 
+ E x E 2 (-2 k ■ p 2 + 3 m 2 - 4 |p x | |p 2 | x x x 2j 
+ 2 |pi| |p 2 | Xx x 2 (k-p 2 + |pi| \p 2 \xx x 2 ))} 
+ C'l [| Pl | 2 | P2 | 2 (-1+x 2 ) (-1+x 2 ) 
x {§E X E 2 - 2 k ■ p 2 + m\ - 6 |pi| |p 2 | x x x 2 \ 

2 (2 E\ E\ + 2k-pxm\-{k- p x f m\ -k-p 2 2 m\- m\ 
2k -px (k-p 2 + mfj |pi| |p 2 | xx x 2 + 2m\ \px\ |p 2 | Xx x 2 



+ El El (-2k ■ p 2 + m 2 e - 6 |pi| |p 2 |xi£ 2 
+ |Pi| |f»2 I xx x 2 (\px\ I F»2 1 xx x 2 (ml - 2 |pi| |p 2 | xx x 2 



+ k ■ p 2 (2m e - 2 |pi| |p 2 | x\x 2 
+ 2 Ex E 2 (k ■ px (k-p 2 + m 2 ) -m A e + k-p 2 (-m 2 + 2 |pi| |p 2 | Xx x 2 
+ |Pi| IP2I ^1 ^2 (-m| + 3|pi| |p 2 | xix 2 ))j } 



(A.2) 
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2. d 2 ^ 0, 
I Fl - 



> l; 

256 7T 



k - pi 



{cf [-3 + 2 E\ E 2 2 + 2{k- Pl ) 2 - Am 2 e k ■ p 2 



Pl 


-4 


IpiI 


IP2I 


P2 


-2 


Pll 


IP2I 



+ |Pi| |f»2 I (-3m 2 xi x 2 + |pi| |p 2 | (l - x\ + (—1 + 3x\ 

+ C' 2 - 

+ Ei E 2 

+ k -pi [—m 2 e +Ak -p 2 - 2 |pi| |p 2 |xix 2 ) 



x 2 



ml + 2 El E 2 + 2(k- Pl ) 2 -3m\k-p 2 
2 1 2k -p! -4|pi| \p 2 \x 1 x 2 



-m„ 



+ 

- 2 

x 

X 



Pil |p 2 | (m%x 1 x 2 + |pi| |p2| 
(k-pi- m 2 e ) k ■ (gi + p 2 ) 

pole— l 
pole+1 



(l-xl+ (-l + 3xj 



E 1 E 2 + k- (pi + p 2 ) + |pi| |p 2 | a;i ^2 - |Pi| |f»2 I yl - 
C' 2 (-ml + k ■ (pi + 2 pa)) + C^ 2 (2m 2 + k ■ ( Pl + 2p 2 ))] } 



x x 2 \[\^x 2 2 



■i -1 



where 



pole 



m 2 e + ExE 2 - k ■ (pi + p 2 ) - |pi| |p 2 | 3:1X2 
~ IPill^l^ 1 -^)^ 1 -^) 



(A.3) 



(A.4) 



3. d 2 ^0, % <1; 



^ = 



256 7T 



{c; 2 [2 £ 2 E% + 2(k- pi) 2 — Am 2 k-p 2 — 3m 



m 2 e — k • pi 

+ £i£ 2 (2/c-pi + 3m 2 -4|pi| |p 2 |xix 2 ) 

+ fc -pi (4/c -p 2 + 4m 2 - 2 |pi| |p 2 | X1X2) 

+ IPil IP2I (-3m e 2 xix 2 + |pi| IP2I (l - A + (-1 + 3xfj a; 2 )) 

+ CJ? [2 J B 1 2 E 2 + 2(A;-p 1 ) 2 -3A;-p 2 m 2 
+ -E 2 (2 fc • pi - m\ - 4 |pi| |p 2 | x x x 2 

+ fc • pi (4 k ■ p 2 - m\ - 2 |pi| |p 2 | xix 2 ) 

+ |Pi| |f»2 I (m\x x x 2 + |pi| |p 2 | (l -x 2 + (-1 + 3x 2 ) x 2 .)) } 



m„ 



(A.5) 



The functions J^ 2 in the same three cases above are obtained from the previous ones 
with the substitutions k — > —k and pi <-> p 2 . 
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B Bremsstrahlung functions 



The functions entering in Eq. (J3.65|) are the following: 



E 
12 



= 512 [C? [{k ■ Pl ) 2 (px ■ p 2 - 2 (ml + k ■ p 2 
- 2m 2 e (m 2 e - k- p 2 +pi-p 2 ) (2 m 2 - k ■ p 2 + p x ■ p 2 
+ 2 k ■ pi (A ml - (k ■ p 2 ) 2 + 3 m 2 e p x ■ p 2 + k ■ p 2 (p 1 -p 2 -2m 2 e 



c' 



2ml -2 ml ((^ ~ V-i) " Pif + (k • Pi) 2 Kra^ — 2 k ■ p 2 +Pi ■ p% 



2 k ■ pi (2 ml + (k ■ p 2 ) 2 - k ■ p 2 (m 2 + p x ■ p 2 ) ) } 
512 [C$ 



(B.l) 



-2 (2m 6 e -Amlk-p 2 + (k- Pl ) 2 (m 2 e + k ■ p 2 ) + m 2 e (k-p 2 ) 2 

+ 2 k ■ pi (3 ml - (2 m 2 e + k ■ p 2 ) (*>2 — Pi • P2)) 

- (Qml - Qm 2 e k ■ p 2 - (k ■ p 2 ) 2 ) P\- p 2 -2m 2 e (p x ■ p 2 ) 2 
+ CJ? [2 TOg - 2 (A; • pi) 2 (ml + k-p 2 j -2k-p 2 (2m A e + k- Pl (k ■ p 2 — m e 
+ Am 2 e (k- p 2 ) 2 + (2 k ■ pi (2 m 2 + k ■ p 2 ) + (& • p 2 ) 2 ) Pi ■ P2 

- 2 m 2 : (pi-p 2 ) 2 ]} , (B.2) 
= 512 [Cl 2 [2mjk-p 2 -m 2 e (k ■ p 2 f - (k ■ Pl ) 2 (m 2 -2 Pl • p 2 ) 

- 2m 4 e pi ■ p 2 + 2m 2 e k ■ p 2 pi ■ p 2 + 2 (k ■ p 2 ) 2 p± ■ p 2 

+ k -pi (k -p2 (3pi ■ p 2 - 4m 2 ) + 2 (m^ + m 2 pi ■ p 2 - 2 (pi • P2) 2 )) 

- 4A;-p 2 (pi -p 2 ) 2 + 2(pi -p 2 ) 3 

+ Ca -mlk ■ p 2 + m 2 e (k ■ p 2 ) 2 + Amlpi ■ p 2 -7 m 2 e k ■ p 2 pi ■ p 2 

+ 2 {k-p 2 f pi -p 2 + {k-pi} 2 (m 2 e + 2pi • p 2 ) + 2 (3m 2 - 2 A; -p 2 ) (pi -p 2 ) 2 

- fc-pi (m^+(7m 2 -3A;-p 2 ) pi-p 2 + 4(pi-p 2 ) 2 )+2(pi-p 2 ) 3 ]} . (B.3) 

The corresponding functions appearing in Eq. (j3.66|) are, instead, obtained from the above 
reported ones by replacing the 4-momentum k with —k. 

After integration over the angle <fi as described in Section 13.2.31 the functions I a for 
the photon absorption and J e0 , I e for the photon emission, appearing in Eq. (|3.72|) . are the 



following^ 



Ei + E 2 ± |k|, £12 = cos^i j2 , s i2 = sin^i^, e 2 = E X E 2 - |pi| |p 2 | x\ x 2 , 
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Ml,2 = #1,2 - |pi,2kl,2 
128 7T Ik 



Mi Ml 



E+ v e 2 /ii /i 2 



a: 2 



3mif4 + m\ n\ (3 m 2 + E+ // 2 ) + m 2 /x 2 /i 2 (3 m 2 - 6 e 2 + 2 £+ /i 2 ) 



+ /i ly u 2 (3m^ + m 2 E+ y u 2 -3e 2 (2m 2 + E+ /i 2 ) 



+ C 



1/2 



c 



'2 



4e 4 /ii/i 2 (/xi +/i 2 ) 



+ (-w 2 Mi + |k| //! /i 2 + /i? /i 2 + /ii /i 2 ) + f4 (|k| /ii - 



-2 m 2 i4 - 2 m 2 + 2 /if /j 2 (2 £+ p /i 2 - m 2 ) + p\ (3 E+ P fx 2 -2m 



+ 2 |px| 2 |p 2 | 2 /ix /i 2 fax + /i 2 ) (1 - x\) (1 - z 2 ,)] } 



(B.4) 



and 



'e(eO) 



'e(eO) 



-64vr k \ /n , 

2-^— (2 7/i +773 

Ml M2 

- 64|k| 

Ml M2 



(ai = 0, a 2 > or a! 7^ 0, a 2 /ai > 1) (B.5) 



[2 (7r — a) 771 — 2 sin a r\ 2 + (tc — a — cos a sin a) 773 



9 sin a + sin 3 a 
6 



r]4 



(ai 7^ 0, — 1 < a 2 /ai < 1) 



(B.6) 



where ai, a 2 , and a have been defined in Section l3\2.3l and we have introduced the following 
functions (£ = for I e while £ = 1 for 7 e0 ): 



Cl 2 



3mJ (/ii + /i 2 ) (/i 2 + /i 2 ) 

+ (M1 + M2) (-4£; 2 E| ) Ui^ 2 +E [ 7 j7 |k|//i/X2(^ + // 2 ) 

+ m 2 (2 e 2 (//i + /i 2 ) (;U 2 - 3 fix fi 2 + // 2 ) + E~ v (/i 2 + /J,) (fij + fiifi 2 + M2)) 

+ C'2 [m H2 (Mi + M2) (-4 -E 2 -E 2 + |k| (/i 2 + /i 2 ) 

+ m\ (2 e 2 (//i + /i 2 ) (/i 2 + /i 2 ) + E~ v (fit + ^ 2 ))' 

+ (Cf + C?) hi M2 [-E~ P e 2 (3 /i 2 + 4 /i! /i 2 + 3 /1 2 ) 

+ 4 (EiE 2 + e 2 ) (/ii +// 2 ) |pi| |p 2 |xix 2 



+ C 



2 e 2 /xi /i 2 - m 2 (//? + //J) 



{^ 2 [-^(e 4 + 3m 2 6 2 + 2 



m. 



+ 3m 2 e \k\(E 1 + E 2 )(fx 1 + l x 2 ) 



+ C' A 



E-- v (e\ E\ - mf) + E~- v (E x E 2 + e 2 ) | Pl | |p 2 | x x x 2 



+ 2 {C' 2 + a 2 ) e 2 \k\(E 1 + E 2 ) (/i! + /i 2 )} 



(B.7) 
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02 



|Pi| IP2M1S2 {6 Cy m 2 /ii n % (/ii +// 2 ) 



+ (Cf + C 



-2 m 2 (/ii + fi 2 ) (fif + i4) + to to (8 e 2 (to + // 2 ) 



+ ^ P (3 y u 2 + 4 y u 1 /i 2 + 3/i 2 ))]} 



+ c 



+ cf 



Ipl 




P2 


SiS 2 




k| 


2 



|-2Cx m\E vD to /i 2 



6 m 2 (Ei + E 2 ) |k| /ii /i 2 (/ii + /i 2 ) + E vV [\2 m 2 e 2 ^1 fi 2 
+ ml (-3/i 2 + /i 2 - 3/x 2 ) 
+ fel 2 + C' 



'/3 



2 £7^ e 2 (3 e 2 /ii // 2 - m 2 e (/i 2 + // 2 ) 
+ 2(E 1 + E 2 ) |k| (/ix + 14) (-4 e 2 Ml M2 + ml (/i 2 + /x 2 ))] } 
= -A(C' 2 + C v 2 ) |pi| 2 |p 2 r / ui / u 2 (/ii + / u 2 )s 2 s 2 



+ c 



Ipl 


2 


P2 


12 2 2 
1 b l b 2 




k| 


2 



{-QC v 2 m 2 e E^Hxto 



+ (Cj? + Cy 2 ) [- 6 E'p e 2 to to + 4 (Ei + £ 2 ) |k| to V2 (fii + /x 2 ) 



C 



2 (d 2 + ^ 2 ) 


£vp 1 


Pi 


3 


P2 


3 Hi H2 s\ s\ 




k 


12 



C Photoproduction functions 

The functions entering in Eq. (J4.77|) are given by 

^ = C v 2 {-32 {k ■ Pl ) 3 m 2 e + Ak- Pl {8mtm 2 + 3m 2 m*) 
+ 4(£;-pi) 2 -8 ml + k ■ p 2 (8 m 2 + 6 m 2 ) + 7 m 2 pi ■ p 2 



+ m 2 e (-7m 2 + I2pi ■ p 2 ) 



m„ 



32 m\ — 16 k ■ p 2 (3 m\ + 2 m 2 m 2 



+ 8 m\ (5 m 2 - 6 p x ■ p 2 ) - pi ■ p 2 + 9 m 2 (m* - 4 m 2 p x • p 2 ) } 

+ Cf {l6 (k ■ pif m 2 e - 8 k - Pl (2 m*m 2 + 3 m 2 m*) 

- 4 (A; ■ pi) 2 (^Ak ■ p 2 m 2 e - Am\ - Q k ■ p 2 m 2 + m 2 e m 2 - 7 m 2 pi ■ p 2 

- m 2 m 2 (16 ml - A (k ■ p 2 - 5 m 2 ) m 2 + 9m*J + m^pi • p 2 } 
+ 2 (Of + C; 2 ) {-8(A;-p 1 ) 2 p 1 -p 2 (fc-p 2 +p 1 -p 2 ) 

+ A(k ■ px) 3 (k ■ p 2 + 2px ■ p 2 ) + 2m 2 (k ■ p 2 + pi ■ p 2 ) [Ak ■ p 2 m 2 e 

+ (4 m 2 + m 2 ) pi • p 2 ) + &; • Pi ^7 (-8 (A; ■ p 2 ) 2 
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+ k -p 2 (5 m 2 - 16 pi -p 2 ) + 4(m 2 - 2p 1 ■ p 2 ) Pi • P2)} , (C- 11 ) 

= Cl 2 {-8k - Pl [4(A;-p 2 ) 2 m 2 + 3 ((k ■ p 2 ) 2 - 2mfj m*-4m 2 e m* 

+ m 2 [32 (k ■ p 2 ) 2 (k ■ p 2 - m\) - 4 (7(k-p 2 ) 2 + 8k-p 2 m 2 e -8m*) m 2 

- 4(3A;-p 2 - 10 m>5 + 9 m°] - [48 m 4 m 2 + 36 m 2 m* + m* 



- 4(A;-p 2 ) 2 (12m 2 + 7m 2 ) 



+ C' 2 {m\ 



Pl -P2 



} 



2\2 2 

y m 7 



-16 (/c • p 2 ) 2 (k ■ p 2 — m 2 ) — 4 (k ■ p 2 — 2 m r 

+ 4 (6 fc • p 2 - 5 m 2 ) m* - 9 m® — 4 k • pi (ml m* 

+ (k ■ p 2 f (-4 m 2 e + 6 m 2 )) - m 2 (-28 (A; • p 2 ) 2 + m*) Pi • p 2 } 

+ 2 (cf + Cf) {fc-pi [4(A;-p 2 ) 3 + 8A;-p 1 (A;-p 2 + m 2 )m 2 

+ hk ■ p 2 mt 1 + 2 [4 k ■ (pi — p 2 ) (k ■ p 2 ) 2 — 8k-pi(k-p 2 + m 2 e ) m 

- ik ■ px + 2k ■ p 2 )mt 1 p x ■ p 2 - 2 [4 (k ■ p 2 ) 2 

- 4(k-p 2 + m 2 e )m 2 -mf f (pi-p 2 ) 2 } , 
= {-4(k-p x ) 2 (8k ■ p 2 m 2 e + (k ■ p 2 - 2m 2 e ) ml) 

+ m 2 \8 (k -p 2 ) 2 ml + 32m\p x ■ p 2 + 5 m 2 (m 2 - 4p x -p 2 )pi ■ p 2 

+ m 2 e (3 m* + 28 m 2 p x ■ p 2 - 48 (pi ■ p 2 ) 2 ) —2k ■ p 2 (4 m 4 

+ 7 m 2 pi ■p 2 + ml(5m 2 + 16p 1 ■ p 2 ))] + 2 k ■ p x \2 (k ■ p 2 ) 2 (8 m 2 + m 2 ) 

- k ■ p 2 (l6m 4 + m 4 + m 2 (22 m 2 - 24 pi ■ p 2 ) + 2m 2 pi • p 2 
+ m 2 (4 m 4 + 7 m 2 pi • p 2 + m 2 (5 m 2 + 16 p 1 ■ p 2 fj } 
+ C'2 {4{k-p l ) 2 k-p 2 {4m 2 e -m 2 )+ml 



(C.12) 



o 2 4 

-3 m e m 7 



+ 4m 4 (m 2 -4pi -p 2 ) + 5m 2 (m 2 - 4pi -p 2 )pi ■ p 2 



- 2 k ■ p 2 (4 m 4 + 7 m 2 pi • p 2 - m 2 (m 2 + 8 pi • p 2 )) 

- 2k -p\ (k -p 2 ) 2 (8m 2 - 2m 2 ) + k ■ p 2 (-8m* + 2m 2 m 2 + m 4 

+ 2 m 2 pi • p 2 ) + m 2 (-4 m 4 - 7 m 2 pi • p 2 + m 2 (m 2 + 8 pi • p 2 )) } 

+ 8 (Cf + C; 2 ) {m 2 p! ■p 2 (k-p 2 +p 1 - p 2 ) (k-p 2 + 2pi • p 2 ) 

+ {k-Pi f ({k -p 2 ) 2 + 2 /c -p 2 pi -p 2 + m 2 pi -p 2 ) 

- k ■ pipi ■ p 2 (2 (k ■ p 2 ) 2 + 2k ■ p 2 pi • p 2 + 3 m 2 pi • P2) } • (C.13) 
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The function J p ^ appearing in Eq. (j4.82j) is obtained after integration over the angle 0. In 
particular, one has 



Iph 

Iph 



128 7r|k| 2 (E 1 -E 2 + E 7 ) 
E 7 (m 2 + 2k- pi) 2 (m 2 -2k- pif 

128 |k| 2 {E l -E 2 + Ey) 
£ 7 (m 2 + 2 k ■ pi) 2 (m 2 -2k- p x f 
9 sin j3 + sin 3/3 



(2 v x + u 3 ) h = 0, 6 2 < 



(C.14) 



[2/91/1+2 sin /5 z/ 2 + (/3 + cos (3 sin /3) z/ 3 



6 



^4 



6i ^ 0, - 1< b 2 /b x < 1 , (C.15) 
and the functions z/j are (xi j2 = cos 6*^2, si >2 = sin^^, e 2 = E x E 2 — |pi| |p 2 | %i x 2 ) 

C' 2 {-gk. pi k-p 2 [2e 6 - Ae 4 k ■ ( Pl - p 2 ) - k ■ ( Pl - p 2 ) ((A; • p x ) 2 + (k ■ p 2 ) 



+ 



X 



+ 



e 2 3(A:-p 1 ) 2 -4A;-p 1 A:-p 2 + 3(fc-p 2 )' 



(k ■ Pl f + e 2 (e 2 -2k- ( Pl - ((* ■ Pl ) 2 + (k ■ p 2 f) + (k ■ p 2 ) L 



+ 16 m 4 e 2 fc ■ pi k ■ p 2 - 8 m 6 e ({k ■ p x ) 2 + (k ■ p 2 ) : 



+ 4 m 



-3 e 4 (A; • pi) 2 + 2 e 6 fc • (pi - p 2 ) + e 2 k • (pi - p 2 ) ((A; • p x ) 2 + (A; • p 2 ) 2 ) 

+ 14 e 4 k ■ pi A; • p 2 — 12 e 2 A; ■ pi /c • p 2 k ■ (pi — p 2 ) 

+ 3 k • pi (fk ■ pi) 2 + (A; ■ p 2 ) 2 ) fc • p 2 - 3 e 4 (fc • p 2 ) 2 - 4 (fc • Pl ) 2 (k ■ p 2 f 

+ 2m 2 ((A;-pi) 3 + e 2 A;-(pi-p 2 )(e 2 -A;-(pi-p 2 )) 

- 2 e 2 ((A; ■ Pl ) 2 + (k ■ p 2 f) + 2 fc ■ Pl k ■ { Pl - p 2 ) k ■ p 2 - (k ■ p 2 ) 3 ) 

- m\ k ■ (pi - p 2 ) (2 e 2 + A; • (p x - p 2 )) - 2 m 6 e k ■ (pi - p 2 ) 

+ 2m 4 [2e 6 + 5e 2 (A;-pi) 2 -8e 4 Av (pi - p 2 ) - 18e 2 k- Pl k-p 2 

+ 4 fc • pi A; • p 2 fc • (pi - p 2 ) + 5 e 2 (k ■ p 2 ) 2 - 2 m\ + m 2 (2 e 2 (e 2 - 2 fc • (p x - p 2 )) 

- 3 ((A-pi) 2 + (A;-p 2 ) 2 ) + 14A;-piA;-p 2 ) -2m 4 (e 2 + 2 fc • (pi - p 2 )] 

- 2 m® 2 e 4 — 3 e 2 k ■ (p x — p 2 ) — k ■ pik ■ p 2 + 5 m 2 A; • (pi — p 2 ) + 2 m 4 

1 8/2 o 2 

+ mile -3 m, 



+ C£ 



12 



2 m 2 + m 2 - 2 e 2 + 2 k ■ ( Pl - p 2 ) 



|2 e 4 (2k-p 2 - m 2 ) (2 A; • p x + m 2 ) + e 2 [-8 k ■ p x k ■ p 2 k ■ {p x - p 2 ) 



Ami 



k ■ pi) 2 + {k ■ p 2 ) 2 ) + 4 fc • pi fc • p 2 ) + 2 m 2 ( (A; • (pi - p 2 )) 



m. 



6 A; • pi k ■ p 2 + 2 m 2 A; • (p x - p 2 )) + 2 m 4 (2 A; • (p x - p 2 ) 
2 ((A;-pi) 2 + (A;-p 2 ) 2 ) + 2m 2 A;-(pi-p 2 )+m 4 ] [2 k ■ Pl (k-p 2 + 
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m,; 



+ m„ 



-2A;-p 2 + 4m 2 + 3m 2 )]} , 



(C.16) 



^2 = |Pi| IP2I s 1 s 2 {C'a {8 A; ■p 1 k 



P2 



4 k ■ pi k ■ p 2 + 6 e 



16 mi 



-8 e 2 k ■ (pi- p 2 ) + 3 ((fc • Pl ) 2 + (A; • p 2 ) 2 ) 
ik- Pl f + e 2 ((A- • Pl ) 2 + (k ■ p 2 ) 2 ) + (A- • p 2 f 
+ k • p\k ■ p 2 k • {p\— p 2 )\ - 16 m 4 A; • pi A; • p 2 + m 2 -4 (6 e 4 A; • (pi - p 2 ) 

6e 2 ((A;. pi ) 2 + (A;-p 2 ) 2 ) + k- Pl k-p 2 (28 e 2 - 13 fc • (p x - p 2 )) - (fc • p 2 f 
+ (k- Pl ) 3 ) + 8m 2 (-2 e 2 fc • (pi - p 2 ) + 3 ((A: • Pl ) 2 + (fc • p 2 ) 2 ) - 2 fc • Pl k ■ p 2 ) 

6 e 4 - 16 e 2 fc • ( Pl - p 2 ) + 5 ((A; • pi) 2 + (k ■ p 2 f) 



+ 8 m 4 A; • (p x -p 2 ) 



-2m 



— 18 k ■ pi k ■ p 2 + 4 m 2 (e 2 — A: • (pi — p 2 )) + 4 rr, 
+ 2m° (4e 2 -3fc-(pi-p 2 )) - m*} 

+ 0? {8k- Pl k-p 2 [6e 4 -8e 2 k-{ Pl -p 2 )+3 ({k ■ Pl ) 2 + {k ■ p 2 )'- 

- 4 A; • pi k ■ p 2 ] - 16 m 2 -(k ■ p{f -2k ■ p x k ■ (pi- p 2 )k ■ p 2 + (k ■ p 2 y 



+ 8m 



3 ((A;.p 1 ) 2 + (A;-p 2 ) 2 ) 



+ e 2 (((A; • Pi) 2 + (k ■ P2 ) 2 ) +6k- Pl k-p 2 )' 

+ 4 k ■ p 1 k ■ p 2 ] + m 2 [-4 ((A; • Pl ) 3 + 6e 4 k-(p 1 - p 2 ) - (k ■ p 2 f 

- 13k-p 1 k-p 2 k-(p 1 -p 2 ) + e 2 (-6(k-p 1 ) 2 + 28k-p 1 k-p 2 -6(k-p 2 ) 2 )) 

+ m\ (32 e 2 k ■ (px - pa) + 12 ((A; • Pi) 2 + (A; • p 2 ) 2 ) + 48 fc • pi fc • p 2 ) 



+ 8 m 4 k ■ (pi - p 2 ) 



-2m 



— m„ 



6 e 4 + 5 (Jfe • pO 2 - 18 k ■ px fc • p 2 + 5 (k ■ p 2 f 
2m 4 -8e 2 (2 k ■ ( Pl - p 2 ) + m 2 )] +2m^ (4 e 2 - 3 k ■ ( Pl - p 2 ) + m 2 ) 

;}} • 

u 3 = 4\p 1 \ 2 \p 2 \ 2 s 2 1 s 2 2 {(c2-2C{ 2 )m 2 e m 2 {2k-(p 1 -p 2 )+m 2 ) 

+ 12 Cy m 2 k-p 1 k-p 2 - (C 12 + C^ 2 ) {4A;-p 1 A:-p 2 (3 e 2 - 2 fc • ( Pl - p 2 

- 2 m 2 ((A; • Pl ) 2 + (A; • p 2 ) 2 ) + m 2 [-6 e 2 k ■ ( Pl - p 2 ) + 3 ((fc • p x ) 2 + (A; • p 2 ) 2 ) 

- Uk-p 1 k-p 2 ] -m 4 (3e 2 -4fc- (pi -p 2 )) , 
1/4 = 4 (CJ? + C^ 2 ) |pi| 3 |p 2 | 3 si s\(2k-p 2 - m 2 ) (2 k ■ px + m 2 ) 



(C.17) 



(C.18) 
(C.19) 
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p 

Figure 8: Self-energy diagram for a photon in a plasma. 

D Electromagnetic excitations in a plasma 

The equation of motion for the 4- vector potential field of a photon in a plasma can be 
written as 

{-k 2 g^ + n M „) A" = , (D.l) 

where is the polarization tensor obtained from the diagram in Figure |H1 (in the Feyn- 
man gauge) jTH] , jH] : 

(D.2) 



with p = p), if = (-E 7 ,k) and E 7 = y|k| 2 + m 2 , -E = y|p| 2 +m 2 . It is useful 
to express n MI/ in terms of form factors by using Lorentz and gauge invariance, namely 
A^IT^i, = PII^ = 0. The most general form of IT^ (assuming parity conservation) results 
to be EH 

IV = n T iv + rizQ^ , (d.3) 

where R^ v and are given in Eqs. (j5.95|) . (|5.94|) . The quantities IT r and TIl may be 
computed from the relations ITl = Q^H^, Ut = R^H^ and the results are 4 : 

4a r | P | 2 (E* g - |k| 2 E 7 / j^+jkp A 
x [F_ (E) + F + {E)} , (D.4) 



4 The expressions in Eas. l|D.4jl and (|D.5|I are obtained by neglecting the term K 4 /A in the denominator 
of Eq. (|D.2(I . As shown in Ref. , this introduces only an error of higher order in a and, moreover, elim- 
inate the effects of the unphysical process 7 — > e + e~, since T1t,l remain real- valued at all temperatures 
and densities as they should be. 
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Ur = 



4a El 



POO 



d\p\ 



P 



71 |k| 2 Jo 
x [F-(E) + F + (E)} , 



E Ulk 



log 



E 1 + u|k| 
E 1 — v\k\ 



E*-\k\ 2 ^ 
E 2 -v 2 \k\ 2 / 



(D.5) 



where v = \p\/E is the electron or positron velocity. The integrals in Eqs. (jD.4|) . ()D.5|) over 
the electron momentum |p| are well approximated ^01 IS] by the following expressions 
that we use in our computations for the neutrino energy loss rates 



2|1,|2 V 



1-G 



k 



+ ^|k| 2 



(D.6) 
(D.7) 



where = uj\/ojp, with the definitions 

u£ = ^ / o °°rf|p| (^-^ 3 ) IPl [F-{E) + F + (E)} 



4a 



— / ^IpI o U 



|p| + 



(D.8) 
(D.9) 



7T JO ' _ ' \3 

can be interpreted as a typical velocity of the electrons in the medium. The function G 
is defined by 



G{x) 



x 



_l 2x 1 — x j / 1 + y/x 



(D.10) 



The dispersion relations -E 7 T,L(|k|) for transverse and longitudinal photon modes are given 
by the locations of the poles in the effective photon propagator which, in the Feynman 
gauge, takes the form 



D 



R 



fJbV 



Q 



fj,v 



(D.ll) 



K 2 -U T K 2 -IL L ' 
The explicit expressions for E lTtL (\k\) are then obtained as the solutions of the implicit 
equations 



E*t,l - |k| 2 = n Tji (E 7TjL , |k|) 



while, near the poles, the scalar parts of the effective propagators are: 



(D.12) 



1 



J T,L 



1 



K 2 - U T r 



2-£' 7 t,l E 1 — E^t,L 



with 



J T,L 



1 + 



T,L 



IE, 



jT,L 



dE, 



(D.13) 



(D.14) 
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By inserting Eqs. (|D.6J) . ()D.7|) into the above equation for Z T ^ Ll we obtain the final ex- 
pressions for the residue functions used in the plasmon decay rate EI] 

2E 2 T (E 2 lT -vl\k\ 2 ) 
ZT = E 2 T (3u£ - 2IL T ) + (E 2 T + \k\ 2 )(E 2 T - vl W) ' (D ' 15) 



Zl = ro a * , ' L • (D.16) 





|k 


I 2 ) 


[3^-(^-«2| 


k| 2 )]n L 



From Ea. (lD~T2j) we can see that £" 7 t > |k| for all |k| while E 1 l > |k| only for |k| < |k| max 
with [IH] 

|k|L, = ^^%I^Qlog(^)-l)[F-(^) + ^)] • (D.17) 

Then longitudinal photon modes can decay into neutrino pairs if their momentum is 
lower than the maximum value |k| max . Note that to the same level of approximation as 



Up 



3 < l .og(i±^)-l 



1/2 

(D.18) 



for Eqs.fjnU, (EH) we hav e [El 

vl \2v* Vl — v 
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